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Abstract. We use the Moderate Resolution Imaging Spec- 
troradiometer (MODIS) satellite aerosol optical thickness 
(AOT) product to assess the impact of reduced swath width 
on global and regional AOT statistics and trends. Along- 
track and across-track sampling strategies are employed, in 
which the full MODIS data set is sub-sampled with var- 
ious narrow-swath (~ 400-800 km) and single pixel width 
(~ 10 km) configurations. Although view-angle artifacts in 
the MODIS AOT retrieval confound direct comparisons be- 
tween averages derived from different sub-samples, careful 
analysis shows that with many portions of the Earth essen- 
tially unobserved, spatial sampling introduces uncertainty in 
the derived seasonal-regional mean AOT. These AOT spatial 
sampling artifacts comprise up to 60 % of the full-swath AOT 
value under moderate aerosol loading, and can be as large 
as 0.1 in some regions under high aerosol loading. Com- 
pared to full-swath observations, narrower swath and sin- 
gle pixel width sampling exhibits a reduced ability to de- 
tect AOT trends with statistical significance. On the other 
hand, estimates of the global, annual mean AOT do not vary 
significantly from the full-swath values as spatial sampling 
is reduced. Aggregation of the MODIS data at coarse grid 
scales (10°) shows consistency in the aerosol trends across 
sampling strategies, with increased statistical confidence, but 
quantitative errors in the derived trends are found even for 
the full-swath data when compared to high spatial resolu- 
tion (0.5°) aggregations. Using results of a model-derived 
aerosol reanalysis, we find consistency in our conclusions 
about a seasonal-regional spatial sampling artifact in AOT. 


Furthermore, the model shows that reduced spatial sampling 
can amount to uncertainty in computed shortwave top-of- 
atmosphere aerosol radiative forcing of 2-3 Wm“^. These 
artifacts are lower bounds, as possibly other unconsidered 
sampling strategies would perform less well. These results 
suggest that future aerosol satellite missions having signif- 
icantly less than full-swath viewing are unlikely to sample 
the true AOT distribution well enough to obtain the statis- 
tics needed to reduce uncertainty in aerosol direct forcing of 
climate. 


1 Introduction 

The direct and indirect effects of aerosols remain the largest 
uncertainties in estimates of the anthropogenic forcing of 
Earth’s climate system (Solomon et al., 2007). Although a 
conceptually simpler problem than the indirect effects of 
aerosols on clouds, the direct effect due to scattering and ab- 
sorption of radiation itself remains poorly constrained owing 
to uncertainty in aerosol loading, temporal and spatial distri- 
bution, and physical properties (Loeb and Su, 2010; Kahn, 
2012). The uncertainty in the anthropogenic direct aerosol 
radiative forcing (DARF) component drives much of the un- 
certainty in overall anthropogenic climate forcing for current 
climate models (Kiehl, 2007). 

Attempts to quantify aerosol properties from satellite ob- 
servations have been made since the 1970s, albeit gener- 
ally with instruments not optimized for observing aerosols. 
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Since the late 1990s, a suite of satellite instruments de- 
signed to measure aerosol properties has helped refine esti- 
mates of aerosol loading, and has contributed some progress 
on retrieving other properties (e.g., absorption, particle size, 
shape, and vertical distribution) (see CCSP, 2009, and ref- 
erences therein). Despite these advances, large uncertainties 
remain, and further reduction of the direct aerosol radia- 
tive forcing uncertainty requires improved satellite coverage, 
as well as integration with in situ observations of aerosol 
type and transport models for synthesis (Diner et ah, 2004; 
Anderson et ah, 2005; Kahn, 2012). 

Given the current state of affairs, design requirements for 
the next-generation, space-based aerosol observing instru- 
ments is a topic of particular importance. For example, in 
2004, the US National Aeronautics and Space Administra- 
tion (NASA), National Oceanic and Atmospheric Admin- 
istration (NOAA), and US Geological Survey (USGS) re- 
quested the US National Research Council (NRC) to “con- 
duct a decadal survey to generate consensus recommenda- 
tions from the Earth and environmental science and ap- 
plications communities” regarding needed space-based and 
ancillary observations to meet future research and opera- 
tional priorities for these agencies (NRC, 2007; hereafter 
referred to as the “decadal survey”). Among the decadal 
survey recommended mission concepts was the Aerosol- 
Cloud-Ecosystems (ACE) mission, which would have as its 
primary goal the reduction of uncertainty about climate forc- 
ing due to aerosols and oceanic CO 2 uptake. A multi-sensor, 
sun-synchronous polar orbiting measurement platform was 
proposed for ACE in the decadal survey. 

In designing a satellite mission, such as the ACE mission, 
there are various trade spaces to consider, an important one of 
which is spatial coverage. Given technological and budgetary 
constraints, trade-offs between the benefit of spatial coverage 
and other desired instrument characteristics must be consid- 
ered, such as the number of spectral and polarized channels, 
relative precision and accuracy, angular and temporal cover- 
age, and pixel size. The challenges of navigating these trade 
spaces for the aerosol problem were recognized in a series of 
papers describing the Progressive Aerosol Retrieval and As- 
similation Global Observing Network (PARAGON, see es- 
pecially Diner et al., 2004; Seinfeld et al., 2004). Anderson 
et al. (2003) addressed the specific question of aerosol spatial 
variability with ground-based and airborne measurements, 
as well as observations from a lidar system flown on the 
NASA Space Shuttle (STS-64), finding important mesoscale 
(40-400 km) structure in aerosol features to be common in 
their data sets. For ACE, the decadal survey explicitly called 
for across-track coverage to resolve this mesoscale structure, 
requirements for which were refined somewhat in NASA’s 
response to the decadal survey ACE concept (Starr et al., 
2010), but those requirements had not been determined from 
any systematic study of aerosol observations, and the afore- 
mentioned conclusions from Anderson et al. (2003) were 
based on spatially and temporally limited data sets. 


Taking advantage of the broader spatial and temporal cov- 
erage afforded by the current space-based, wide-swath Mod- 
erate Resolution Imaging Spectroradiometer (MODIS, see 
Sect. 2.2), Geogdzhayev et al. (2013) investigated the impli- 
cations of along-track only sampling on statistics of global 
monthly mean aerosol optical thickness (AOT). AOT is rele- 
vant to the question of DARF, as it is a proxy for aerosol col- 
umn loading, and because to first order it determines DARF. 
Hansen et al. (1995) suggest that a change in the global mean 
AOT of 0.01 corresponds to a climatically important change 
in the global mean radiative forcing of 0.25 Wm“^. This 
can be compared with the 0.5 zb 0.4 W m“^ Intergovernmen- 
tal Panel on Climate Change (IPCC)-stated uncertainty in the 
magnitude of the anthropogenic DARF component (Solomon 
et al., 2007). Other analyses suggest that the actual uncer- 
tainty is far larger than the IPCC estimate (McComiskey et 
al., 2008; Loeb and Su, 2010). Arguing from their analysis 
of MODIS AOT retrievals, Geogdzhayev et al. (2013) con- 
cluded that narrow along-track-only sampling, only one pixel 
wide, is suflhcient to constrain the global, monthly mean aver- 
age AOT, and is therefore adequate for climate applications. 
They did not, however, address implications of reduced spa- 
tial sampling on regional AOT, or the uncertainty in decadal- 
scale AOT trends. 

In this paper we re-examine the question of the impact of 
spatial sampling for a sun-synchronous polar orbiting satel- 
lite on the statistics of AOT. We approach this question 
using two complementary data sets. The first data set ex- 
plored is the AOT retrieval product from the MODIS instru- 
ment aboard the NASA Aqua spacecraft, as in Geogdzhayev 
et al. (2013). Several different along-track and across-track 
sampling strategies are applied to the MODIS observations 
in order to sub-sample the data set. The second data set is 
derived from a newly available aerosol reanalysis performed 
in the NASA Goddard Earth Observing System (GEOS-5) 
model. GEOS-5 is an Earth system model that includes an 
online aerosol module, and is run in our case using assim- 
ilated meteorology and invoking assimilation of AOT from 
satellite observations. In Sect. 2 we describe and present re- 
sults of our sampling study as applied to the MODIS data set. 
In Sect. 3 we investigate the question of spatial sampling in 
the context of the GEOS-5 model results. Section 4 presents 
a discussion of our results and our conclusions. 


2 Spatial sampling of the MODIS data set 

2.1 A conceptual illustration of the spatial sampling 
problem 

At any given time, nature presents us with a particular three- 
dimensional spatial distribution of clouds and aerosols, as 
well as the attendant variability in particle microphysical 
characteristics, surface reflectivity, and solar illumination. 
The passive satellite instrument retrieval problem amounts 
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to inverting a meaningful geophysical quantity (e.g., AOT) 
from this complexity, given a limited set of measured param- 
eters (e.g., backscattered spectral reflectance). Our hypothe- 
sis is that the ability to tease out the climatically significant 
portion of this signal for synoptically important events de- 
pends in part on the spatial and temporal coverage of the ob- 
serving system. In this paper we focus on spatial coverage as 
determined by the sensor’s viewing-swath width. 

We illustrate the spatial coverage aspects of the problem 
conceptually in Fig. 1 . The global aerosol field is represented 
with a familiar painting (Fig. la). In a perfect observing sys- 
tem we would see the true aerosol field at all times and lo- 
cations, but lacking that we are constrained by the realities 
of our instrument’s spatial and temporal coverage, as well 
as by the challenges inherent in the retrieval problem itself. 
We show this using a single day’s coverage of the aerosol 
field retrieved by the MODIS instrument aboard the Aqua 
spacecraft. Superimposed on the painting are the orbital gaps, 
clouds, and bright surfaces where the MODIS “dark target” 
aerosol retrieval is not applied. Note the large white area in 
the shape of the Saharan and Arabian deserts over the fig- 
ure’s left eye (Fig. lb). For purposes of our illustration, we 
are making the assumption that our “true” scene is time in- 
variant, as it would take MODIS 24 h to achieve the cov- 
erage shown in Fig. lb. If we narrow our view further, as- 
suming something like the spatial coverage of the NASA 
Multi-angle Imaging Spectroradiometer (MISR) instrument, 
we obtain the scene in Fig. Ic, and further reduction of the 
coverage to a single pixel wide “curtain” along the MODIS 
track yields Fig. Id. This sampling construction is formally 
developed in Sect. 2.3. Figure 1 illustrates that very different 
pictures of the “true” scene emerge depending on the spatial 
coverage of the observing system. In what follows, we quan- 
tify the impact of spatial coverage characterizing the time 
varying global and regional field of AOT. 

2.2 The MODIS data set 

We use aerosol observations from the space-based MODIS 
instrument for this part of our study. MODIS provides near- 
global, daily AOT retrievals over land and ocean surfaces. 
There are two MODIS instruments, both in sun-synchronous 
polar orbits. MODIS on the Terra satellite has been opera- 
tional since early 2000 and has a daytime equator crossing 
time of about 10:30 a.m. local time (LT) at the center of its 
swath. MODIS on the Aqua satellite has been operational 
since mid-2002 and has a daytime equator crossing time of 
about 01:30 p.m. LT. At the nominal orbit altitude of 704 km, 
the MODIS instruments observe a swath about 2300 km wide 
along their ground tracks. The MODIS orbit is such that the 
ground coverage is repeated exactly every 16 days. AOT is 
retrieved in the daytime portion of the MODIS orbit under 
cloud-free and glint-free conditions using separate aerosol 
retrieval algorithms for ocean (Tanre et al., 1996, 1997) and 
land (Kaufman et al., 1997; Levy et al., 2007a, b). 
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In our analysis, we use the land and ocean AOT retrievals 
from the MODIS Aqua instrument, valid at 550 nm, from the 
Collection 5 MODIS algorithm products (Remer et al., 2005, 
2008; Levy et al., 2010). The retrievals are made at a nom- 
inal 10 km X 10 km spatial resolution at nadir. A quality as- 
surance (QA) flag is reported for each retrieval, indicating 
its estimated level of confidence as a valid result, from tests 
performed during the retrieval process. QA flags range from 
0 (lowest confidence) to 3 (highest confidence). In order to 
retain the highest quality MODIS data, in what follows we 
use only the highest confidence (QA = 3) retrievals over land, 
and require QA>0 over ocean (Remer et al., 2008; Levy et 
al., 2010; Breon et al., 201 1). The uncertainty in the MODIS 
AOT (r) product is characterized such that one standard devi- 
ation (66 %) of the retrievals fall within Ar = ±0.03 zb 0.05 r 
over the ocean and Ar = ±0.05±0.15r over land relative 
to the AOT from coincident ground-based AERONET sun 
photometer network observations (Remer et al., 2005). The 
AOT data set thusly obtained is representative of the clear- 
sky (cloud-free) aerosol distribution, and it should be under- 
stood that a portion of our subsequent spatial sub-sampling 
(see next section) is entwined with the spatial sampling of 
clouds. 

2.3 Methodology for along-track sub-sampling of AOT 
from the MODIS full swath 

Our spatial sampling strategies are illustrated in Fig. 2, 
which shows an example over-ocean scene comprising a sin- 
gle MODIS Aqua swath. Along-track spatial sub-sampling 
(along the MODIS orbital track) is illustrated in Fig. 2a. We 
consider the AOT retrieved across the MODIS full swath 
(FS), as well as several sub-sampled swaths in which we re- 
tain only the relevant portions of the full swath. Four nar- 
row swaths are chosen to approximate the "^380 km wide 
swath of the Multi-angle Imaging Spectroradiometer (MISR, 
on the Terra spacecraft. Diner et al., 1998), labeled: N1 - 
east side of full swath, N2 - adjacent to Nl, N3 - adjacent 
to N2 and in center of swath, and N4 - west side of swath. 
We also consider a “mid-width” swath (MW) with coverage 
between the narrow and full swath composed of the union of 
Nl andN2. To approximate the “curtain”-like sampling of an 
instrument such as the Cloud Aerosol Lidar with Orthogo- 
nal Polarization (CALIOP, aboard the CALIPSO spacecraft. 
Winker et al., 2010), we consider the samplings Cl, C2, C3, 
and C4, which are extracted at the center of the Nl, N2, 
N3, and N4 swaths, respectively. As will be explained be- 
low, the MODIS glint artifact particularly impacts the N3 
and C3 samplings. We emphasize that in what follows, we 
are using only MODIS AOT retrievals, sub-sampling the full 
data set along the indicated narrow and curtain swaths. The 
sampling strategies are summarized in Table 1. Figure 2 also 
shows across-track sampling and sampling of global model 
data that will be described and used in Sects. 2.5 and 3.2, 
respectively. 
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Figure 1. Conceptual illustration of the spatial sampling problem. Nature presents us with a “true” seene (a). The truth is sampled aeeording 
to “full-swath” sampling of the MODIS instrument on the Aqua spaeeeraft (b), a “narrow” sampling (c), and a single pixel wide “eurtain” 
sampling (d). For purposes of this illustration we are reeovering only parts of the “true” image that had valid aerosol retrievals on 5 June 2010 
from the MODIS over-oeean and “dark target” land retrievals. 


The individual retrievals after sampling are aggregated 
onto regular latitude-longitude spatial grids typical of the 
grids used in global aerosol transport models. We will in 
this paper consider two grid resolutions: 0.5° x 0.625° and 
10° X 10°. For each, the grid-averaged AOT is 

n 

E ■ (li 

%id = ;; . (1) 

II li 
! = 1 

where tt are the 1 through n individual AOT retrievals falling 
into the grid box and qi is the QA value assigned to each re- 
trieval. Our aggregation is thus QA weighted. Over land we 
have only retained QA = 3 retrievals, based on the MODIS 
Aerosol Product Data Quality Statement. The aggregation 
is performed daily. The temporally averaged (e.g., monthly, 
seasonal, annual) AOT at a grid box is 

m 

^ %id,; • nj 

’ ( 2 ) 

7=1 

where m is the number of days to average, rgUd,; is the grid 
average value at day j from Eq. (1), and nj is the number of 
retrievals used to make tgud,;. This aggregation and weight- 
ing strategy is the same as in Remer et al. (2008) and Colarco 
et al. (2010). 


2.4 Along-track sub-sampling results 
2.4.1 The sub-sampled AOT 

The sub-sampled MODIS Aqua AOT data are analyzed for 
the years 2003-2012. Figure 3a shows, for example, the year 
2010 annually averaged AOT from the full-swath MODIS 
Aqua retrievals over both land and ocean using the aggrega- 
tion strategy given by Eqs. (1) and (2) and the 0.5° x 0.625° 
aggregation grid. The spatial patterns of the main aerosol 
features are evident: the Saharan dust and Asian pollution 
and dust outflow plumes, the biomass burning activity over 
southern Africa and South America, the pollution plume over 
China, the band of high AOT in the southern ocean, and a 
region of high AOT over western Russia where a signifi- 
cant biomass burning anomaly occurred in 2010 (Witte et al., 
2011 ). 

For comparison, the corresponding 2010 annual mean 
AOT of one of our narrow (Nl) and two of our curtain (Cl 
and C3) samplings are shown in Fig. 3b, c, and d, respec- 
tively. As seen in Fig. 2, Nl and Cl are on the eastern edge 
of the MODIS swath. C3 is down the center of the swath 
and is particularly impacted by the sunglint artifact, and so 
has relatively poor retrieval sampling over the tropical ocean. 
Many of the features apparent in the full-swath annual mean 
in Fig. 3a are still apparent in the sub-samples: the biomass 
burning plumes over South America and southern Africa, the 
Asian outflow across the northern Pacific, Saharan dust trans- 
port across the North Atlantic, and dust and anthropogenic 
pollution over India and China. On the other hand, the shapes 
and apparent magnitudes of these features are clearly differ- 
ent, and certain features are notably absent, particularly the 
Saharan dust plume in the C3 sampling (mostly in the glint 
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Figure 2. Spatial sampling strategies used in this paper illustrated for an example portion of an over-oeean orbit traek of the Aqua spaeeeraft 
on 5 June 2010. (a) Along-traek sampling illustration. The eolored dots indieate the loeations of the MODIS AOT retrievals, with the 
grey dots indieating the full MODIS swath (MO). Overlaid on the grey dots are different eolors for our various along-traek sampling 
strategies (N1 = light blue, N2 = orange, N3 = magenta, N4 = light green. Cl = dark blue, C2 = dark red, C3 = deep purple, C4 = dark green, 
and MW = eombined N1 and N2 swath). The light-grey shaded areas on the left and right side of the figure are outside the swath, while 
the eentral white region (labeled “glint”) is where no aerosol retrievals are made due to glint. Remaining patehy white areas are where 
aerosol retrievals were not made due to elouds. (b) As in (a) exeept illustrating aeross-traek sampling following Geogdzhayev et al. (2013). 
(c) Along-traek sampling strategy for MERRAero model output. Dark grey shading is the MODIS Aqua full swath. Light and dark blue are 
respeetively the narrow and eurtain width sampling along the Aqua ground traek. Orange and dark red are respeetively the narrow and eurtain 
width sampling along the CALIPSO ground traek. Note the absenee of eloud or glint features in the model sampling. 


Table 1. Summary of spatial sampling strategies illustrated in Fig. 2 and summary of temporal averaging approaehes. 


Sample name 

Sample width 

Full swath (FS) 

~ 2300 km 

Mid-width (MW) 

~ 800 km 

Narrow (4 variants: Nl, N2, N3, N4) 

~ 380 km 

Curtain (4 variants: Cl, C2, C3, C4) 

~ 10 km (width of MODIS pixel) 

Averaging strategy 

Proeedure 

Sample-then-average 

Per orbit, sample the MODIS full swath at the indieated sub-swath 
Aggregate sub-sample to spatial grid 
Average aggregates to the desired time period 
(e.g., monthly, seasonal, annual) 

Average-then-mask 

Per orbit, aggregate the MODIS full swath to spatial grid 
Average to the desired time period 
Use “sample-then-average” result for relevant 
sub-sample/temporal average to retain or 
exelude grid boxes visited in sub-sample 


region) and the high AOT features over the southwest United 
States in the Cl and N1 samplings. Other features, such as 
the 2010 Russian fires (Witte et al., 2011), still exhibit sim- 
ilar magnitude of AOT across sampling strategies, but their 


spatial coherence is quite different. Grey shading indicates 
places where no retrievals are available at all during the year. 

We emphasize that in Fig. 3 the approach is “sample- 
then-average,” and so is done on a “per-orbit” basis (see 


www.atmos-meas-tech.net/7/2313/2014/ 


Atmos. Meas. Tech., 7, 2313-2335, 2014 



2318 


P. R. Colarco et al.: Impact of satellite viewing-swath width 





AOT [550 nm] (Annual mean year 2010) 
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Figure 3. MODIS Aqua year 2010 annual mean AOT using the sampling and aggregation strategy in Eqs. (1) and (2) shown for various 
along-traek spatial samplings: (a) full swath, (b) Nl, (c) Cl, and (d) C3. The grey shading indieates loeations where no MODIS AOT 
retrievals were made during the year. 


Table 1). Only the MODIS retrievals that could have been 
sampled are pulled from the full-swath data set, then aggre- 
gated, and then finally averaged. This “sample-then-average” 
approach is how time averages are typically calculated from 
polar orbiting satellite data sets. We make this point to dis- 
tinguish from a different sampling approach discussed later 
(Sect. 2.4.2). 

Figure 4 shows the years 2003-2012 time series of global, 
annual mean AOT over both land and ocean for each of our 
sampling strategies generated with a similar procedure to 
what is shown in Fig. 3 (grey lines). Also shown in Fig. 4 are 
the global annual mean AOT resulting from other sampling 
approaches, to be discussed in Sects. 2.4.2, 2.5, and 3.1. 

The full-swath annual mean AOT varies between about 
0.13 and 0.14 over the ocean and between about 0.16 and 
0.18 over the land, similar to the multi-year analysis pre- 
sented in Remer et al. (2008). We compare the global, an- 
nual mean AOT of our various sampling strategies to the 
full-swath AOT. Over ocean, except for the N4 and C4 sam- 
plings, the global, annual mean AOT is within 0.01 of the 
full-swath value. Over land, most of the sampling strategies 
differ from the full swath by more than 0.01 at some point in 
the time series, with Nl and Cl notably underestimating the 


global, annual mean AOT relative to the full swath. The ex- 
cursions of the individual sampling strategies global, annual 
mean AOT from the full-swath mean can be as much as 15 % 
of the full-swath magnitude. This is important, because if the 
narrow-swath or curtain-like sampling cannot reproduce the 
basic statistics of the full-swath AOT even at the global and 
annual scales, the question of whether we can rely on this 
measurement strategy for narrowing the uncertainties in key 
aerosol properties and their impacts on climate must be as- 
sessed quantitatively. 

There is, however, a significant caveat to our along-traek, 
“sample-then-average” results as presented in Fig. 4. Al- 
though the differences between any given sub-sample and the 
full-swath AOT certainly contain a component related to the 
spatial sampling, errors and uncertainties in the MODIS re- 
trievals themselves also contribute to the differences shown. 
In particular, the MODIS AOT retrieval has a sensor view- 
angle dependency (Levy et al., 2010). That is, if the aerosol 
loading is homogeneous across the MODIS swath, different 
AOT values will under some circumstances nevertheless be 
retrieved in different positions across the swath, owing to this 
angular artifact. However, the characteristic of that artifact 
as a function of view angle, sun angle, position on Earth, 
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Figure 4. A time series for the years 2003-2012 of the global, annual mean MODIS Aqua AOT over oeean (a) and land (b). Grey shaded 
lines are the MODIS along-traek sample-then-average approaeh. Red lines are the MODIS along-traek average-then-mask approaeh. Blue 
lines are the aeross-traek sampling following Geogzhayev et al. (2013). Blaek lines are the MERRAero along-traek sampling. Different lines 
within eaeh eolor group are the various sampling eandidates explored. The bottom panels (c, d) in eaeh is the differenee of the individual 
sub-sampled averages from the respeetive full-swath average. Individual MERRAero sub-samples are differeneed relative to the MERRAero 
full-swath (Aqua) sampling. Others are relative to the MODIS Aqua full-swath sampling. On top panels only, green dashed line is the data 
assimilation (DA) grade MODIS AOT assimilated in MERRAero (see text). 


surface reflectance, etc., is not well understood. In an earlier 
study on the sampling question posed here (Colarco et al., 
2012), we attempted to correct for this dependency by exam- 
ining a data set of MODIS-AERONET collocations sorted by 
view geometry, similar to what is shown in Levy et al. (2010) 
(see, for example, their Fig. 10). This proved challenging. 
The collocation data set was relatively small and was only 
available where AERONET sites are located. The latter point 
made it difficult to evaluate the view-angle dependency of 
the MODIS AOT retrievals, especially over ocean. The data 
set would have been smaller still for determining these view- 
angle dependencies on a seasonal or regional basis. For these 
reasons, we could not separate view angle from spatial sam- 
pling differences when the full-swath and sub-sample AOT 
data sets were compared using the sample-then-average ap- 
proach so far employed, so from this point forward we take a 
different approach to evaluate the impact of swath width on 
global AOT statistics. 

2.4.2 Application of an observability criterion 
to the MODIS data set 

Although we cannot correct the MODIS observations for the 
view-angle dependency with confidence, we can investigate 
the question of observability: what are the characteristics of 


the observations that are not made in a given sub-sample, and 
how does this impact the derived AOT statistics? 

Figure 5 complements Fig. 3. It shows the 2010 full- 
swath annual mean AOT from Fig. 3a, but only in grid boxes 
where the indicated sub-sampling strategy had no valid an- 
nual mean AOT (i.e., in grid boxes where the sub-sample 
either never visited because of coverage limitations, or else 
never encountered a good AOT retrieval because of algorith- 
mic issues when it did overfly the grid box). We show this for 
Nl, N3, Cl, and C3 sub-samples. This is revealing. The N1 
and N3 figures show that, like the full-swath sampling, the 
narrow-swath sampling permits retrieval of AOT over most 
points on Earth at least once during the year. However, in 
some places where the full-swath sampling would have made 
relatively few observations, the narrow-swath sampling pro- 
vides no observations at all. (See Fig. 6, which shows the 
number of aerosol retrievals in each grid cell considered for 
several of the sampling strategies discussed in this paper.) 
These regions are generally where seasonal changes in sur- 
face brightness due to vegetation (e.g., the US southwest, the 
Sahel) or seasonal snow cover (the Tibetan plateau) make 
retrieval difficult and thus less frequent. On the other hand. 
Fig. 5 illustrates something qualitatively different when only 
curtain-like sampling is obtained: it is clear that much of the 
Earth is never visited at all under this sampling. 
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Figure 5. MODIS Aqua full swath year 2010 annual mean AOT shown only at points never sampled by the indieated sub-sample swath: 
(a) Nl, (b) N3, (c) Cl, and (d) C3. 


In our analysis of observability, to reduce the issue of the 
view-angle artifact discussed above we create what we call 
our “average-then-mask” strategy (Table 1). First, we con- 
struct monthly, seasonal, or annual mean maps of the AOT 
from the full-swath data, effectively sampling the location at 
all viewing geometries obtained by MODIS. Second, we cre- 
ate masks that exclude the grid boxes never observed by each 
particular sampling strategy over the relevant averaging pe- 
riod. Finally, we apply the masks to the aggregated maps of 
the full-swath AOT. This “average-then-mask” strategy is in 
contrast to the “sample-then-average” strategy described in 
Sect. 2.3. 

The results of this method provide a view of the features 
each sampling strategy can observe, and estimates of the 
mean AOT differences that are unbiased by scan angle ar- 
tifacts. But it also represents a much richer data set than 
could be obtained from an instrument having similar retrieval 
capabilities to the full-swath MODIS but having only nar- 
row or curtain sampling. As a result, this method reduces 
significantly the difference in AOT variability measured by 
the different sampling strategies compared to the difference 
obtained using the “sample-then-average” method. This re- 
duction in the variability is illustrated also in Fig. 4 for the 
“average-then-mask” approach (red lines). It is clear that in 


the global, annual mean, the AOT for the sub-samples gener- 
ally differs from the full-swath value by much less than 0.01 . 
However, the “average-then-mask” approach for the curtain- 
like cases provides spatial sampling that is much greater than 
would ever be acquired by an actual curtain instrument, be- 
cause the aggregated samples come from different parts of 
the broad MODIS swath. So although this approach mini- 
mizes the view-angle bias, it includes far more of the broad- 
swath data than would be available from a curtain instrument. 

2.4.3 An illustration of the average-then-mask sampling 
approach 

Using the “average-then-mask” method, we emphasize the 
issue of observability further in Fig. 7, where we zoom in on 
the key aerosol features in the region surrounding the trop- 
ical Atlantic Ocean, focusing on the seasonal AOT for the 
period July- August- September 2010. We show the MODIS 
full-swath seasonal mean AOT at grid cells both where Cl 
and C3 do and do not observe (i.e., the union of the two Cl 
sampling images. Fig. 7a and c, yields the full-swath sea- 
sonal mean for this region, as does the union of the two C3 
sampling images. Fig. 7b and d). Here the spatial gaps in 
the curtain-like sampling become more apparent, and visual 
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Figure 6. Number of MODIS Aqua AOT retrievals made per 0.5° x 0.625° grid box for the entire year 2010 as used to eompose the (a) 
full-swath annual mean shown in Fig. 3a and for the (b) LI, (c) Nl, and (d) Cl sub-samples. Light grey shading indieates loeations where 
no retrievals were available. 


inspection of the figures reveals differences in the patterns 
of the aerosol features seen. For example, the AOT features 
over the Nile River valley are nearly absent in the Cl and C3 
sampling (Fig. 7a and b) but are readily apparent when look- 
ing at the grid cells unobserved by these sampling strategies 
(Fig. 7c and d). Likewise, the biomass burning plume over 
South America is more structured, and has a much greater 
range of values, for the grid cells unobserved by Cl and C3 
(again. Fig. 7c and d). Most readily apparent is the wide 
equatorial belt over the ocean, encompassing the Saharan 
dust plume, where the C3 sampling is almost completely ab- 
sent due to glint (Fig. 7b versus d). Even for the Cl sampling, 
where the ocean glint is not an issue in this case, the South 
African biomass burning plume is also missing some of the 
highest- AOT regions when the observed and unobserved grid 
cells are compared (Fig. 7a and c). 

A similar analysis is presented in Fig. 8 over the Asian re- 
gion for the March- April-May 2010 seasonal average. Here 
we show only the Cl sub-sample masking. Aerosol features 
over Iraq, Iran, Turkmenistan, Afghanistan, northern China, 
and the Sichuan Basin in central China are almost completely 
unobserved by the Cl sampling (Fig. 8a), the seasonal AOT 
maxima and minima are lost, and the pattern of the main 
Asian outflow over the northern Pacific is less well defined. 


2.4.4 The AAOT spatial sampling artifact in 
seasonal-regional analysis 

For AOT trend and regional climate impact studies, quan- 
titative differences matter. We assess the quantitative dif- 
ferences produced by different sampling strategies for sev- 
eral regions exhibiting major aerosol features as highlighted 
with white boxes in Figs. 7 and 8. We introduce a diagnostic 
metric, AAOT, which we call the AOT spatial sampling ar- 
tifact. For a given spatial and temporal average, AAOT is 
the spread in the spatiotemporal mean AOT among all of 
the sampling strategies considered. So, for each of the re- 
gions highlighted in Figs. 7 and 8, we compute the “average- 
then-mask” seasonal-regional mean AOT for our time se- 
ries for all of the along-track sampling strategies consid- 
ered. AAOT is then the difference between the maximum 
and minimum value among the sampling strategies for each 
seasonal-regional point: 

AAOT = Tmax '^min- (3) 

Because the glint significantly impacts the sampling in the 
C3 and N3 sub-samples for certain regions, we exclude those 
sub-samples from the AAOT calculation in all cases (i.e., 
AAOT is determined from the remaining eight sampling 
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Figure 7. Full-swath seasonal (July- August- September 2010) MODIS Aqua AOT over the tropieal Atlantie Oeean. The full-swath seasonal 
mean is masked to show only grid eells where the Cl and C3 sub-samples do (a, b) and do not (c, d) have a seasonal mean value. Figure 7a 
and b illustrate the “average-then-mask” seasonal mean AOT. The white boxes outline regions diseussed in the text. 
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Figure 8. As in Fig. 7, but for the Cl sub-sampling mask for Mareh-April-May 2010 over Asia. The full-swath seasonal mean AOT is shown 
both where the Cl sub-sample does (a) and does not (b) have a valid seasonal mean. The white boxes outline regions diseussed in the text. 


strategies). This restriction proves especially important for 
the southern Africa, African dust, Nile River, Southeast Asia, 
and Asian outflow regions, where glint would otherwise con- 
found the results. 

Figures 9 and 10 present, respectively, the A AOT spatial 
sampling artifacts for South America and the Indo-Gangetic 
Plain, two of the eight regions highlighted in Figs. 7 and 8 


(the other regions, discussed here, are shown in the Supple- 
mentary material). We focus here on the along-track sam- 
pling artifact (red lines), and will discuss the results from 
other sampling approaches later in Sects. 2.5 and 3.2. To 
highlight the differences between curtain-like and narrow- 
swath sampling, we also show the AAOT for the full swath. 
Cl, C2, and C4 samplings only (dotted red line) and for 
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South America 
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Figure 9. Seasonal-regional sampling artifaets AAOT and 
Aforeing for South Ameriea. Red lines are the MODIS along-traek 
AAOT sampling artifaet using the average-then-mask approaeh 
(solid: all samples exeept C3 and N3; dotted: full swath, Cl, C2, 
and C4 samples only; dashed: full swath, Nl, N2, and N4 sam- 
ples only). Blue line is the MODIS aeross-traek AAOT sampling 
artifaet (solid for the “sample-then-average” approaeh, dashed for 
the “average-then-mask” approaeh). Blaek line is the along-traek 
AAOT sampling artifaet from the MERRAero results. Grey shading 
(eorresponding to right-hand vertieal seale) is the Aforeing artifaet 
from MERRAero results. Also shown are the MODIS full-swath 
mean AOT (bottom, solid line), the range of the along-traek sub- 
sample AOT values about the full-swath mean (bottom, shading), 
and AAOT as a fraetion of the full-swath AOT (bottom, dashed). 
Finally, we show the eorrelation eoeffieients between AAOT and 
the full-swath seasonal-regional mean AOT, as well as the AAOT 
as a fraetion of the full- swath mean. 

the full swath, Nl, N2, and N4 samplings only (dashed red 
line). For all, we additionally show the full-swath AOT value, 
the range of the sub-sample AOT values around the full- 
swath value, and the magnitude of AAOT as a fraction of 
the seasonal-regional full-swath AOT. Finally, the corre- 
lation coefficient of the AAOT with the full-swath AOT and 
the fraction of the full swath are also indicated. 

We refer to AAOT as the “sampling artifact”, as it shows 
the uncertainty in the seasonal-regional AOT due to limited 
spatial sampling, relative to the full swath. We note that for 
all regions the AAOT sampling artifact is highest for the 
curtain-like sampling (dotted red line), and so drives the sam- 
pling artifact when all sampling strategies are included (solid 
red line). We also emphasize that AAOT as computed here 
is actually a lower bound on the possible AOT spatial sam- 
pling artifact. If additional along-traek sampling strategies 
were considered they could possibly contribute to increasing 
the spread between the maximum and minimum seasonal- 
regional AOT, but they would not reduce the spread over 
what already is computed. 


Indo-Gangetic Plain 



Figure 10. As in Fig. 9, but for the Indo-Gangetic Plain. 


The South America region (Fig. 9) shows significant an- 
nual and inter-annual variability in the full-swath AOT, with a 
peak AOT of between 0.2 and 0.4 typically occurring in July- 
August- September (JAS) or October-November-December 
(OND) associated with seasonal biomass burning. This peak 
is modestly correlated (r^ = 0.25) with the AAOT, which 
can be as high as 0.06. Interestingly, AAOT is uncorrelated 
with its fractional comparison to the full-swath AOT (^laction 
in bottom panel of Fig. 9, which is the correlation of bottom 
panel dashed red line and top panel solid red); although, as 
a fraction of the full-swath AOT the AAOT typically peaks 
at 40 % and can be as high as 60 %. Thus, for South Amer- 
ica, the uncertainty in AOT due to sampling may be as much 
as 0.06, comprising ~ 15 % of a base magnitude as high as 
about 0.4, and can also represent uncertainties as great as 
60 % in the regional AOT when AOT is lower. 

Southern Africa (Fig. Sla in the Supplement) is another 
region affected by seasonal biomass burning, with peak AOT 
of about 0.4 occurring in JAS. The peak AAOT is at most 
0.03 and is weakly correlated with the full-swath AOT (r^ = 
0.14), but much more strongly with the fractional contribu- 
tion (r^ = 0.72). 

For African dust (Fig. Sib in the Supplement) the AAOT 
is small (approximately 0.01) and is consistently less than 
about 5% of the magnitude of the full-swath seasonal- 
regional AOT. In other words, for the African dust region, 
the average-then-mask sampling does not significantly im- 
pact these AOT statistics. For the Nile River (Fig. Sic in the 
Supplement) the AAOT is at most about 0.05 and is modestly 
correlated {r^ — 031) with the full-swath seasonal-regional 
mean AOT signal. The full-swath mean AOT has a seasonal 
signal, varying between about 0.2 and 0.4 in magnitude, and 
the sampling artifact may be as much as about 20 % of the 
foil-swath value. 

Turning to Asia, for the Indo-Gangetic Plain (Fig. 10), the 
peak values of AAOT are as high as 0.1 but are uncorrelated 
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with the full-swath AOT, which itself peaks in magnitude at 
about 0.5. The sampling artifact may thus be as much as 
about 30% of the full-swath signal. In China (Fig. Sid in 
the Supplement) the A AOT is as high as 0.09 and is some- 
times as large 20 % of the full-swath mean AOT, which itself 
varies between about 0.3 and 0.6 in magnitude. In the South- 
east Asia (Fig. Sle in the Supplement) and Asian outflow 
(Fig. Slf in the Supplement) regions the peak A AOT val- 
ues are 0.05 and 0.015, respectively. For Southeast Asia, this 
sampling artifact can be as large as 20 %, but is mostly less 
than 10% of the full-swath signal. The contribution to the 
Asian outflow signal is negligible, with sampling introduc- 
ing an uncertainty of only about 5 % at most for a full-swath 
AOT that peaks above 0.4 in magnitude. 

In summary, with the “average-then-mask” approach, dif- 
ferences presented are due solely to sampling, as we are 
only comparing the gridded full-swath seasonal mean val- 
ues with sub-samples of itself, and cross-swath anomalies 
are implicitly removed by the averaging. This “average-then- 
mask” approach incorporates much greater sampling than ac- 
tual reduced-swath instruments can obtain - about three-to- 
four times more samples for the narrow-swath, and about 
16 times more samples for the curtain. Nevertheless, sig- 
nificant qualitative and quantitative differences still appear 
in the seasonal-regional average AOT distributions; minima 
and maxima do not capture the extreme values, and some re- 
gional features are entirely missed. Due to the much greater 
sampling included in the “average-then-mask” data, results 
presented in Figs. 4, 9, 10, and SI in the Supplement are 
significantly more favorable than would be produced for real 
instruments having such spatial sampling characteristics, and 
thus the sampling artifacts presented in this section are effec- 
tively lower bounds. The overall magnitude of the sampling 
artifact is largest for the curtain-like sub-samples, as might 
be expected. The nature of this artifact is such that in some 
regions (South America, Indo-Gangetic Plain, China) it can 
be as large as 20-60 % of the full-swath AOT signal or as 
great as 0.1 in AOT magnitude. 

2.4.5 Along-track sub-sampling and trends in aerosol 
optical thickness 

In the previous section, we showed artifacts that are intro- 
duced in the seasonal-regional mean AOT when the full- 
swath data are sub-sampled. Here we investigate how spatial 
sampling affects the ability to detect statistically significant 
AOT trends. 

Our approach follows the trend analysis presented in 
Zhang and Reid (2010), which employs the statistical tools 
of Weatherhead et al. (1998) to assess confidence levels in 
the derived trends. Briefly, a linear model is fit to the monthly 
mean AOT time series at a grid box. A first-order autoregres- 
sive “noise” model characterizes the residual of the observed 
time series from the linear model. The slope of the linear fit, 
(D, is the trend in the time series, and the standard deviation 


of the trend, is defined in terms of the variance of the 
residual noise model (see Eq. 2 in Weatherhead et al., 1998). 
Where the ratio >2, the trend is statistically signifi- 

cant at the 95 % confidence level (Weatherhead et al., 1998). 

AOT trends and statistical significance maps are produced 
for our 10-year time series of MODIS AOT data. As in Zhang 
and Reid (2010), the AOT monthly mean time series is de- 
seasonalized before the linear model is fit, because the sea- 
sonal aerosol signal is so large in many parts of the world. 
We considered the trends both in terms of the monthly means 
composed from our “average-then-mask” and “sample-then- 
average” methodologies. In Figs. 11 and 12, we show the 
“sample-then-average” trends and statistical significance, re- 
spectively, for the MODIS full swath and several of our sub- 
samples. The “sample-then-average” method more realisti- 
cally represents the data that would be acquired by a narrow- 
swath or curtain instrument, but is subject to the uncertainty 
associated with the MODIS view-angle artifact, and so may 
contain scan angle biases in the AOT field itself that could 
alias the magnitude of the derived AOT at some locations. 
However, this will not affect the statistical significance of the 
derived trends as long as whatever scan angle artifacts exist 
do not vary over time for a given sub-sample of the MODIS 
swath. The high calibration stability of the MODIS instru- 
ments (Xiong et al., 2006) supports this assumption, although 
a calibration drift in certain MODIS channels does affect the 
Collection 5 MODIS AOT data that has since been corrected 
for Collection 6 (Levy et al., 2010, 2013). For the purpose 
of the current study, we are concerned primarily with differ- 
ences in the statistical significance of the trends that can be 
derived for various distributions of samples. 

Figure 1 1 shows the AOT trend for several of our sampling 
approaches. The full swath (Fig. 11a) shows a strong negative 
trend in AOT over the Amazonian region in South America, 
in Siberia north of Mongolia, and across Alaska, and mod- 
erate negative trends across the eastern United States and 
Canada and the western North Atlantic Ocean, Europe and 
the Mediterranean, in the Gulf of Guinea and off the west 
coast of northern Africa, and in the western Pacific around 
the Maritime Continent. Strong positive trends are apparent 
in the Arabian Sea, across India, and in the Bay of Bengal, in 
Iraq, off the western coast of southern Africa, across Sudan 
and Ethiopia, near Beijing in eastern China, in eastern cen- 
tral Argentina, and in eastern Siberia and across the northern 
Pacific Ocean. Moderate positive AOT trends are seen in the 
western United States and Canada, over southern Africa, and 
more generally across northern Asia. Except as noted previ- 
ously, the oceans generally have no trend in AOT or else a 
weakly positive trend. The locations and signs of our com- 
puted trends are generally similar to Zhang and Reid (2010, 
their Fig. 7a). Our trends differ from theirs primarily in the 
Pacific west of Mexico, where we show a slight positive trend 
and they show a weak negative trend, and as well in the west- 
ern Pacific and near the Maritime Continent where they find 
a weak positive trend and we find a weak negative trend. We 
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Figure 11. Trend for the 10-year (2003-2012) time series of MODIS Aqua AOT. We show the trend for the (a) full swath, (b) Nl, (c) Cl, and 
(d) LI samplings. Grey areas are loeations with either no valid retrievals or where the time series has fewer than 12 months and months _i 
pairs. 




note that Zhang and Reid (2010) restricted their analysis to 
over-ocean regions only, derived trends from a different time 
series (2000-2009), and used their “assimilation-grade” ver- 
sion of the MODIS AOT product, which is quality controlled 
as described in Zhang and Reid (2006), and so is a somewhat 
different data set than considered here. 

The Nl sampling trends (Fig. 1 lb) are also similar in pat- 
tern and sign to the full-swath trends, but differences from the 
full swath are more clearly visible, including a stronger posi- 
tive trend associated with the southern African biomass burn- 
ing plume and a more strongly negative trend across central 
western Africa and in northeastern Asia. The lesser coverage 
associated with the Cl sampling makes the trends harder to 
discern for that case (Fig. 1 Ic), although the overall patterns 
of increasing and decreasing trends are again fairly consistent 
with the full swath. The other narrow and curtain-like sam- 
ples have similar trend patterns and magnitudes (not shown), 
but differ in detail, and the N3 and C3 samples have poor 
coverage over the tropical oceans. We also considered the 
mid-width (MW) sampling trend (not shown), finding it to 
be similar in behavior to the full-swath trends. (The across- 
track LI sampling trend will be discussed in Sect. 2.5.) 


The differences in the trend magnitudes between our sam- 
pling approaches are not unexpected. Zhang and Reid (2010) 
found, for example, weaker magnitude for trends from MISR 
observations than for MODIS, and attributed this difference 
at least in part, if not entirely, to the lower spatial coverage 
of MISR. Additionally, our “sample-then-average” approach 
can affect the magnitudes of the trends, due to MODIS view- 
angle biases discussed previously. Our focus is thus on our 
ability to assign statistical significance to whatever trends ap- 
pear in the maps. 

In Fig. 12 we present the spatial distribution of statisti- 
cal significance for the trends shown in Fig. 11. For the full 
swath (Fig. 12a) our analysis shows that the computed trends 
are significant at the 95 % level broadly across the tropical 
Pacific Ocean, in the Arabian Sea and Bay of Bengal, in the 
Mediterranean, and then across Sudan and Ethiopia and into 
the western Indian Ocean. Our patterns for significance are 
again very similar to those of Zhang and Reid (2010) (their 
Fig. 7b), except that their broad region of significance be- 
tween southern Africa and South America is much less pro- 
nounced in our analysis. Over land we find statistical sig- 
nificance in the full swath for southern India, near Beijing, 
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Figure 12. As in Fig. 11, but showing the statistical significance for the trends shown in Fig. 11. Regions colored red show regions with 
statistically significant trends at the 95 % confidence level. 


across the central United States, in Argentina, and across por- 
tions of the biomass burning region in Amazonia. 

For the MW sampling (not shown) the over-ocean patterns 
of significance are nearly identical to the full swath, but over 
land there are notable differences, with MW indicating no 
significance in the derived trends over Amazonia, in China, 
or in the central United States. The regions of significance in 
the trends over India and in the Sudan and Ethiopia region are 
also much reduced in area. This reduction in areal extent of 
significance patterns worsens for the N1 sampling (Fig. 12b), 
with significance essentially gone over Sudan and Ethiopia, 
and as well being much reduced over Argentina. The pat- 
terns over ocean are still generally similar to the full swath, 
but the individual regions are less coherent. For the Cl sam- 
pling (Fig. 12c), the statistical significance at the 95 % confi- 
dence level is essentially gone, with nothing identifiable over 
Sudan, Ethiopia, India, or Argentina, sparse identification of 
significance over other land area (especially at high northern 
latitudes), and only a hint of significance in the tropical Pa- 
cific and in a few other ocean regions. The patterns of signif- 
icance for the other narrow and curtain-like samplings (not 
shown) are similar to the N1 and Cl shown in Fig. 12, re- 
spectively, though different in detail. N2 and C2 have some- 
what better coverage over the oceans. N3 and C3 - again. 


because of the glint - show poor coverage over the oceans. 
C3 in particular has far worse coverage over the ocean than 
Cl. Although the reduction in spatial sampling reduces the 
number of grid cells in which statistical significance can be 
identified versus the full swath, it is possible that at least in 
broad regions - for sufficient averaging areas - we can nev- 
ertheless obtain statistically significant trends even for the 
curtain sampling. We address this quantitatively in the next 
sub-section. 

For comparison, the trends and statistical significance 
were also computed from monthly means constructed using 
our “average-then-mask” method. Recall that this provides 
far greater sampling than the actual along track sampling 
would provide. The patterns, signs, and magnitudes of the 
trends (Fig. S2 in the Supplement) are quite similar to those 
shown in Fig. 11. The effect of the more favorable “average- 
then-mask” approach is most notable for the patterns of sta- 
tistical significance (Fig. S3 in the Supplement), which are 
broadened considerably for the N1 and Cl samples in this 
approach. Still, considerable spatial sampling gaps remain, 
especially over India, the Arabian Sea, Sudan and Ethiopia, 
Brazil, and the central US, where statistical confidence in the 
derived trends is not assured in the narrower samplings. 
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2.5 Across-track sub-sampling results 

The recent paper by Geogdzhayev et al. (2013) is of par- 
ticular relevance to this study, as they provided a similarly 
motivated analysis of the MODIS AOT data. Their approach 
was to develop sub-samples by aggregating individual scans 
across the MODIS swath. They argued that averaging across 
the MODIS swath removed the view-angle artifact when 
compared to the full set of MODIS observations, versus a 
comparison to along-track sampling (i.e., samples similar to 
our Cl-4 sub-samples). 

We implemented this sampling approach in the same 
framework as the along-track samplings discussed earlier, 
selecting five evenly spaced across-track sub-samples (LI, 
L2, L3, L4, and L5, with the “L” standing for “latitudi- 
nal”, see Fig. 2b). The year 2010 annual mean AOT for 
the LI sub-sample is shown in Fig. 13. When compared 
with the full-swath annual mean AOT (Fig. 3a), we see a 
lot of “noise” (small-scale variability) in the AOT field for 
the LI sub-sample. Consistent with the earlier discussion of 
our along-track sub-samples, there are important aerosol fea- 
tures missed by this sampling, including the South American 
biomass burning plume and the Russian fires. Nevertheless, 
when the global, annual mean AOT is compared to the full- 
swath AOT, there is essentially no difference between any 
of the latitudinal sub-samples and the full swath (Fig. 4, blue 
lines). The annual means shown in Fig. 4 are for the “sample- 
then-average” methodology applied to this data set. Our re- 
sult is consistent with Geogdzhayev et al. (2013), who also 
used the “sample-then-average” methodology. Even though 
this approach was shown not to work so well for along- 
track sampling (grey lines in Fig. 4), we find that if any- 
thing there is less variability among the sub-samples for the 
across-track sub-sampling than for the along-track with the 
“average-then-mask” approach applied (i.e., spread in blue 
lines is smaller than the spread in the red lines in Fig. 4). 
When the “average-then-mask” approach is applied to the 
across-track samplings the spread in the global, annual mean 
AOT collapses even further. 

When considering the seasonal-regional statistics, how- 
ever, it is clear a spatial sampling artifact still remains in the 
across-track sampling, as might be expected from the small- 
scale variability in the map of global AOT (Fig. 13). Fig- 
ures 9, 10, and SI in the Supplement includes the across- 
track AAOT spatial sampling artifact, which can be di- 
rectly compared to the along-track results discussed earlier. 
Both the “sample-then-average” and “average-then-mask” 
variants of this are shown for the across-track sampling. 

In general, the “average-then-mask” across-track sampling 
artifact is similar in magnitude to the along-track values, 
in some regions on average less and in others on aver- 
age greater. In all regions, however, there are compara- 
ble magnitude spikes in this sampling artifact. By contrast, 
the “sample-then-average” across-track sampling artifact is 
typically much larger than the along-track or across-track 
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Figure 13. As in Fig. 3d, the year 2010 annual mean AOT gridded at 
0.625° X 0.5° horizontal resolution for the LI aeross-traek MODIS 
sampling. 


“average-then-mask” values, corresponding to its relatively 
poorer spatial coverage. 

In Figs. I Id and 1 2d we show, respectively, the AOT trend 
and statistical significance pattern for the LI across-track 
sub-sample. The global distribution of the sign of the trends 
is generally consistent with the full-swath data set (Fig. 1 la), 
but there are considerable differences in coverage. The full- 
swath observations have hundreds-to-thousands of observa- 
tions per year informing a given grid box (Fig. 6a), whereas 
the LI -type sampling has at most a few dozen (Fig. 6b). The 
relatively poor coverage for the LI sampling at this resolu- 
tion renders the trend statistically insignificant almost every- 
where (Fig. 12d). This is also true for the other latitudinal 
sub-samples (not shown). The particular areas of coverage 
and trend magnitudes differ somewhat among the different 
latitudinal sub-samples, but in all cases there is almost no 
ability to assign statistical significance. The result is the same 
when the “average-then-mask” approach is applied to calcu- 
lating the trend (Fig. S2d in the Supplement) and significance 
(Fig. S3d in the Supplement). 

For completeness, we performed this same trend analysis 
at a coarser 10° x 10° spatial aggregation, compatible with 
the resolution of the analysis performed in Geogdzhayev et 
al. (2013). Here, the sampling of the MODIS data is per- 
formed on the coarse resolution grid, monthly means are 
formed, and the trends are calculated. The AOT trends and 
the map of 95 % statistical significance for the full swath, 
LI, Nl, and Cl samplings are presented in Figs. 14 and 15, 
respectively. Results may be compared with Figs. 11 and 
12. The trend magnitude and sign at the coarser aggrega- 
tion resolution are qualitatively similar to those derived at the 
higher aggregation resolution, and at the coarse resolution 
the trends across the sub-samples are generally similar to one 
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Figure 14. Aerosol trends for the (a) full swath, (b) Nl, (c) Cl, and (d) LI samplings at 10° x 10° aggregation resolution. 


another. The assignment of significance to a detected trend is 
of course statistically more robust at the coarser spatial res- 
olution, since relatively more of these larger grid boxes have 
valid monthly means at the coarser spatial resolution. Thus, 
unlike what was seen at higher spatial resolution (Fig. 12c 
and d, respectively), at 10° x 10° spatial resolution it is pos- 
sible to assign statistical significance more broadly for the Cl 
and LI samples (Fig. 15c and d). Using the more favorable 
“average-then-mask” approach, the coarse resolution aggre- 
gated trends (Fig. S4 in the Supplement) and statistical sig- 
nificance (Fig. S5 in the Supplement) are even more alike 
across the sub-samples. 

3 Spatial sampling of a global aerosol transport model 

In analyzing spatial sampling issues in the context of the 
MODIS AOT data set, as was done in Sect. 2, instrumental 
and algorithmic artifacts in the data set contribute uncertainty 
to our conclusions. To that end, we seek to also apply our 
analysis machinery to a data set that is as independent as pos- 
sible from the MODIS product. It seems natural, then, to in- 
vestigate these issues in the context of a global aerosol trans- 
port model, where - whatever its shortcomings - at least it is 
not impacted by, e.g., view-angle artifacts, calibration drift, 
sun glint, or even - necessarily - clouds. However, we must 


also consider that an aerosol transport model presents issues 
of its own. There is a considerable literature within the Ae- 
roCom community, for example, that documents the perfor- 
mance and biases with current aerosol transport models (e.g., 
Textor et al., 2006; Kinne et al., 2006; Schulz et al., 2006; 
Huneeus et al., 2011; Myhre et al., 2013; Stier et al., 2013; 
Randles et al., 2013). To the issues explicitly called out in 
those references, we add that the ability to represent aerosol 
spatial and temporal variability on scales smaller than 100s 
of kilometers and months in time is particularly relevant to 
the question being addressed here. Climate models with grid 
boxes several 100s of kilometers in size seem to us not to be 
the appropriate tool. Fortunately, we have access to a recent 
aerosol reanalysis produced in a relatively high spatial res- 
olution global aerosol transport model. Although we cannot 
here resolve all the possible issues in applying this aerosol 
reanalysis to our question, or for that matter provide a com- 
prehensive validation of the reanalysis, we hope that we offer 
sufficient confidence to proceed with our central question. 

3.1 The GEOS-5 MERRA aerosol reanalysis 
(MERRAero) 

The NASA Goddard Earth Observing System Model 
(GEOS-5) is an Earth system model maintained by the 
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Figure 15. As in Fig. 14, but for the 95 % statistical significance interval. 


NASA Global Modeling and Assimilation Office (Rienecker 
et al., 2008). GEOS-5 contains components for atmospheric 
and oceanic circulation, land surface processes, ocean bio- 
geochemistry, radiation, and atmospheric composition and 
chemistry. An aerosol module derived from the Goddard 
Chemistry, Aerosol, Radiation, and Transport (GOCART, 
Chin et al., 2002) model has been incorporated in the GEOS 
system (Colarco et al., 2010). The GOCART module treats 
the life cycles - including sources, sinks, and chemistry 
- of tropospheric sulfate, carbonaceous, dust, and sea salt 
aerosols as a series of non-interacting tracers (see the Sup- 
plement). 

In addition to functioning as an Earth system model, a pri- 
mary focus of GEOS-5 itself is as a tool to provide reanal- 
yses to NASA instrument teams and the broader scientific 
community. The present atmospheric reanalysis performed in 
GEOS-5 is the Modem-Era Retrospective Analysis for Re- 
search and Applications (MERRA, Rienecker et al., 2011), 
which provides a global, three-dimensional (0.5° x 0.666° 
latitude by longitude, 72 vertical levels from the surface 
to approximately 85 km), six-hourly atmospheric reanaly- 
sis for the period 1979-present. Augmenting MERRA is 
our own aerosol reanalysis, so-called MERRAero, covering 
the period 2002-present. MERRAero is a “replay” of the 
GEOS-5 system driven by the MERRA atmospheric fields. 


with the GOCART aerosol module mn online within the 
system, and including assimilation of AOT derived from 
MODIS Aqua and Terra (Kessner et al., 2013; Buchard et 
al., 2014; Sessions et al., 2014). The AOT data assimilated 
in MERRAero are not the MODIS Collection 5 aerosol 
retrievals themselves, but rather result from a neural net- 
work regression of the Collection 5 reflectances to inde- 
pendent AOT measurements from the Aerosol Robotic Net- 
work (AERONET, Holben et al., 1998), a ground-based, 
global network of sun photometers (see the Supplement). 
The global, annual mean of this derived “data assimilation “ 
(DA) grade MODIS Aqua AOT that is assimilated in MER- 
RAero is shown in Fig. 4 (dashed green line), and is quite 
obviously lower in magnitude than the MODIS Collection 
5 AOT (e.g., red and blue lines), in part due to the more 
conservative cloud clearing being applied. For comparison, 
Zhang and Reid (2010) also show a lower magnitude for their 
DA grade version of the MODIS AOT data set than is found 
for the operational MODIS retrieval algorithms. 

MERRAero is run globally at 0.5° x 0.625° latitude by 
longitude horizontal spatial resolution, and we analyze the 
period 2003-2012, to cover the same analysis period ap- 
plied to the MODIS data. We look at the AOT simulated 
in MERRAero, and although this AOT has an imprint of 
MODIS information on it due to the aerosol assimilation. 


www.atmos-meas-tech.net/7/2313/2014/ 


Atmos. Meas. Tech., 7, 2313-2335, 2014 





2330 


P. R. Colarco et al.: Impact of satellite viewing-swath width 


we draw from the model the simulated AOT immediately 
prior to the aerosol assimilation step, which means that ef- 
fectively we are looking at model AOT fields that have not 
seen the MODIS information for approximately 24 h. Over- 
all, MERRAero has an AOT bias with respect to AERONET 
of —0.022, an RMS error of 0.148, and an correlation co- 
efiicient of 0.629. These statistics are improved over versions 
of the model that do not invoke aerosol assimilation. Further 
description and evaluation of the MERRAero system and the 
aerosol assimilation methodology are a work in progress (see 
also the Supplement). 

3.2 Results of sampling analysis 

The gridded nature of the model fields lends itself to a some- 
what different approach to sub-sampling than was developed 
for the MODIS retrievals. We use an orbital sampling exten- 
sion to the OpenGrADS software package (http://opengrads. 
org/doc/udxt/orb/orb.html) which exploits the readily avail- 
able two-line element (TEE) orbital parameters from a wide 
array of operational spacecraft in order to develop an orbital 
masking tool. This orbital mask is “fiown” through the grid- 
ded model output, and a number of different sub-samples are 
accumulated. Figure 2c illustrates how this is applied to sam- 
pling MERRAero fields for a single orbit along the Aqua 
track (a MODIS-like full swath, a MISR-like narrow swath, 
and a CALIPSO-like “curtain” track), as well as for an orbit 
of the approximately parallel nearby Aqua spacecraft (a nar- 
row swath and a “curtain” track). Note that the operational 
TEE is not used in this package, but rather a canonical TEE 
for the spacecraft, and so the orbital sampling is not precisely 
the same as what the Aqua (or Aura) spacecraft actually fiew. 
Additionally, we note that in contrast with Fig. 2a and b, there 
are no cloud or glint features in Fig. 2c. MERRAero is sub- 
sampled without regard to these issues as the model is ca- 
pable of computing an AOT in all grid boxes regardless of 
cloud cover or other conditions that affect satellite viewing. 
The AOT from MERRAero is the all-sky AOT; that is, it is 
derived from the grid box mean aerosol and humidity fields. 

The global, annual mean AOT time series from the MER- 
RAero sub-samples are shown in Fig. 4 (black lines). MER- 
RAero over-ocean AOT tends to be lower than the MODIS- 
derived sub-samples, and there is a trend toward decreasing 
AOT apparent in MERRAero that is not seen in the MODIS 
data. Furthermore, over most of the time series MERRAero 
has a higher over-ocean AOT than the DA-grade MODIS- 
derived AOT being assimilated. For now, we attribute these 
aspects of MERRAero to the fact that, despite having an 
assimilated, satellite-derived AOT product, MERRAero re- 
mains an aerosol transport model. So, for example, as data 
sets are not yet available to drive anthropogenic aerosol emis- 
sions in the model for the full time period being simulated, 
there is an imprint of model error present that the data as- 
similation evidently does not entirely correct. Over land the 
MERRAero global, annual mean AOT is more comparable in 


magnitude to the MODIS data, but the inter-annual variabil- 
ity is quite different. Again we emphasize that MERRAero 
is not subject to surface boundary conditions affecting the 
MODIS retrievals, and the over-land MERRAero averages 
include AOT values over bright desert and snow-covered sur- 
faces. For both land and ocean the MERRAero global, annual 
mean AOT for all sub-samples are within 0.01 of the MER- 
RAero full-swath averages. 

The MERRAero AAOT spatial sampling artifact is shown 
in Figs. 9, 10, and SI in the Supplement. It is generally the 
case that the MERRAero sampling artifact is smaller than 
what we derived from MODIS. This difference, however, is 
not always great. Seasonal peaks in the sampling artifact in 
South America, for example, can be as large in MERRAero 
as in MODIS, and are even somewhat larger in MERRAero 
for the African dust region. On the other hand, in China the 
seasonal peaks in sampling artifact are at most about half 
the magnitude of what was found in the MODIS analysis. 
In MERRAero, the sampling artifacts are small for the Asian 
outfiow and Southern Africa regions, similar to what MODIS 
showed, whereas the Nile River region notably has a large 
excursion above the magnitude found in MODIS in April- 
May-June 2007. This excursion is traced to a May 2007 dust 
event in the Sahara that produces a high regional AOT for the 
curtain sampling and a low regional AOT for the full-swath 
sampling, highlighting the variability of small sources and 
the potential impact on seasonal-regional AOT statistics. A 
similar feature is seen in Southeast Asia for the JEM season 
in 2010, in this case attributable to a biomass burning feature 
preferentially sampled by the curtain sampling. 

Finally, we exploit the fact that we are analyzing results 
of an aerosol transport model here. In the model run cycle, 
a radiative transfer code is called every time step, and top 
of atmosphere (TOA) and surface shortwave and long-wave 
radiative fiuxes are computed. For each grid column the radi- 
ation call is made several times, both including and removing 
the effects of aerosols, and similarly the effects of clouds. 

In Figs. 9, 10, and SI in the Supplement, we also show the 
sampling artifact as computed in terms of the TOA all-sky 
shortwave fiux with and without aerosols (this is the TOA 
shortwave aerosol direct radiative effect), where the magni- 
tude of the sampling artifact Aforcing is again the difference 
in the seasonal-regional mean between the samples having 
the highest and lowest mean values, analogous to Eq. (3). 
As we are sampling only during the daytime portion of the 
satellite orbit, we approximate the diumally averaged radia- 
tive forcing sampling artifact by dividing the actual artifact 
computed from the instantaneous values by a factor of four. 
We note that the all-sky TOA aerosol radiative forcing is neg- 
ative (cooling) over all regions considered except the Nile 
River region, where the forcing is positive owing to the ex- 
ceedingly bright desert surface. The Aforcing artifact is by 
construction defined positive. 

In all regions except the Asian outflow region (where 
the AOT sampling artifact is small), the Aforcing sampling 
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artifact has seasonal peaks in excess of 1 W m“^, and it ap- 
proaches 3Wm“^ in the Nile River region. There is not 
always a strong or obvious correlation between the MER- 
RAero Aforcing and AAOT sampling artifacts, reflecting 
the complexity of the forcing being a convolution of aerosol 
loading with the distributions of clouds and the varying, un- 
derlying surface properties. For example, the same Aforcing 
sampling artifact can be computed under clear-sky conditions 
(i.e., removing the effects of clouds), in which case the forc- 
ing itself is more strongly correlated with the AOT in all re- 
gions. 


4 Discussion and conclusions 

We have used MODIS data and results of a global aerosol re- 
analysis to investigate the impact of spatial sampling on the 
statistics of AOT. From MODIS data, we showed significant 
differences in the global, annual mean AOT derived from 
along-track sampling (Fig. 4). The “sample-then-average” 
approach employed initially, however, could not disentangle 
the spatial sampling artifacts (which we are most interested 
in isolating) from the sensitivity of the MODIS AOT retrieval 
to viewing geometry. 

Subsequently, we considered the observability problem: 
where the sub-sample could have obtained aerosol retrievals, 
where it could not, and where - compared to the full- 
swath values - differences in the regional and seasonal AOT 
are found. The observability-based “average-then-mask” ap- 
proach (Sect. 2.4.2) mitigates biases associated with loca- 
tion in the MODIS swath, but greatly increases the sam- 
pling compared to an actual instrument having a narrower 
swath, because the full-swath MODIS instrument obtains 
much more frequent observations of any given location than 
an actual narrow-swath instrument would. This approach 
yielded global, annual mean AOT values that were insignif- 
icantly different from the full-swath AOT values, in contrast 
to the originally employed “sample-then-average” method; 
nevertheless, important regional differences remained. The 
AAOT spatial sampling artifact was typically small for our 
more ocean-influenced regions (e.g., African dust, Asian out- 
flow), but could be as large as 0.1 in the seasonal AOT aver- 
age for high-loading, near-source regions such as China and 
the Indo-Gangetic Plain. As a percentage of the full-swath 
seasonal-regional mean AOT, the sampling artifact was as 
large as 60% (South America), and was in many places 
around 20 % (China, Indo-Gangetic Plain, Nile River). In al- 
most all cases, the magnitude of the sampling artifact was 
largest for the curtain-like sampling, with smaller artifacts 
inferred when the narrow-swath sampling was compared to 
the full swath, as might be expected. 

Further analysis of the MODIS data using the across-track 
sampling method suggested in Geogdzhayev et al. (2013) 
likewise yielded global, annual mean AOT values insignif- 
icantly different from the full-swath values, but again with 


significant regional differences, often exceeding the sam- 
pling artifact for our along-track sampling. 

Finally, the MERRAero reanalysis results sampled along 
hypothetical orbital tracks also yielded global, annual mean 
AOT values insignificantly different from the full-swath sam- 
pling of the MERRAero AOT field. The MERRAero AAOT 
spatial sampling artifacts were generally similar in magni- 
tude to the along-track MODIS sampling artifact. From the 
MERRAero results, we also computed a Aforcing sampling 
artifact as the range in the seasonal-regional mean all-sky 
TOA aerosol forcing across all sampling strategies. In most 
regions this Aforcing artifact had seasonal peaks larger than 
1 W m“^, and in several regions it had seasonal mean peaks 
as large as 2-3 W m“^. The Aforcing artifact is distinct from 
the AAOT artifact in that it bears an imprint of surface and 
cloud variability in addition to aerosol loading variability, 
and so reveals a different aspect of the spatial sampling is- 
sue explored here, i.e., the ability to constrain the climati- 
cally critical aerosol direct radiative forcing (Solomon et al., 
2007). 

Our analysis was focused on seasonal-regional spatial 
sampling artifacts, which capture aspects of the mean aerosol 
distributions. Folded into these mean field properties is vari- 
ability occurring on much smaller spatiotemporal scales, in- 
cluding extreme events (i.e., volcanic eruptions, wildfires) 
that could be missed depending on the spatial sampling. 
The “average-then-mask” method employed for our MODIS 
along-track results would carry events to the sub-samples 
more completely than would be the case for the “sample- 
then-average” approach that more realistically represent the 
sampling of an actual instrument. Because of the aforemen- 
tioned limitations of the “sample-then-average” data set, we 
do not fully assess here the sub-regional distribution of AOT 
and how that is impacted by spatial sampling, but this is an 
issue that needs further consideration. 

We additionally investigated our ability to detect sta- 
tistically significant trends in aerosol features as a func- 
tion of spatial sampling. We find that the full-swath trends 
in our study are similar to the “contextually less-biased” 
assimilation-grade trends in Zhang and Reid (2010), sug- 
gesting that although contextual bias can be an issue, it 
probably does not diminish the applicability of our conclu- 
sions. Although the signs of the trends were similar for the 
various sampling strategies employed, magnitudes were in 
some places quite different. This is attributable in part to the 
MODIS view-angle bias, but also to differences in the spa- 
tial coverage. Again, most places on Earth are simply never 
observed with curtain-like sampling, including some major 
aerosol source regions. When the MODIS data set is grid- 
ded at the highest spatial resolution used here, the reduced 
spatial coverage has a profound impact on the ability to as- 
sign statistical significance to the trends (Fig. 12). For ex- 
ample, even the narrow N 1 sampling could not assign sig- 
nificance at the 95 % confidence level to any decadal-scale 
trends over Amazonia or the central United States, and had 
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reduced confidence in western Africa and India. The pat- 
terns of significance were even less coherent for the curtain- 
like sampling. Analyzing the trends with the “average-then- 
mask” approach yields improvement in the ability to assign 
statistical significance in the narrower samples, but there are 
still broad regions where statistical significance cannot be as- 
sured without the full-swath sampling. Although the narrow- 
swath sampling fares better than curtain sampling, without 
the context of a full swath imager’s observations, there is lit- 
tle confidence in even these derived trends, a conclusion sim- 
ilar to one obtained by Zhang and Reid (2010). The conclu- 
sion is similar for the across-track sampling of the MODIS 
data, where - again - although the derived trend is similar 
in magnitude to the full-swath-derived trend, the ability to 
assign statistical significance is compromised by the poorer 
spatial sampling. On the other hand, aggregation of the sub- 
sampled AOT fields to coarser spatial resolution (10° x 10°, 
Figs. 14 and 15) drew out many common features in both 
the trend magnitude and statistical significance across sam- 
pling strategies, suggesting that narrow and even single-pixel 
wide sampling are sufiicient for qualitative trend detection at 
these scales. Note that for sufficiently coarse spatial and tem- 
poral aggregation, all global sampling strategies should con- 
verge to the same values, regardless of sampling frequency. 
However, for finite time series, the ability to capture regional 
effects and to detect trends depends heavily on sampling fre- 
quency. 

We note that the MODIS data set does not capture all as- 
pects of the actual aerosol field, in part due to contextual 
limitations of the measurement technique, such as the lack 
of diurnal observations and the inability to retrieve AOT un- 
der and in the immediate vicinity of clouds (e.g., Zhang and 
Reid, 2009). Others have identified contextual issues with the 
MODIS data set over land (Hyer et al., 201 1) and in coastal 
sites (Anderson et al., 2013). Other measurement approaches 
would enhance retrieval of aerosol properties over bright land 
surfaces and address these issues. Future aerosol instruments 
will undoubtedly improve upon MODIS data set in these and 
other respects, such as providing enhanced information about 
aerosol single scattering albedo and particle size, as well as 
other important drivers DARF (e.g., Loeb and Su, 2010). 
It was not our intention here to demonstrate the benefit of 
those enhanced capabilities, but rather to investigate the lim- 
its imposed on the measured AOT by one aspect of any fu- 
ture measurement strategy, its spatial coverage. That our re- 
sults are broadly consistent between the analysis of MODIS 
data - however imperfect - and an aerosol transport model - 
subject to its own biases, but biases which are different than 
MODIS ’s - increases confidence in the results. 

Our main conclusions are thus that 

- Spatial sampling is not a driver of uncertainty in the 
global, annual mean AOT, for the sampling strategies 
considered. This finding is robust for the “average-then- 
mask” along-track and across-track MODIS sampling. 


as well as for the MERRAero sampling, but not so for 
the “sample-then-average” approach, which aside from 
MODIS swath-related artifacts, is a better representa- 
tion of the sampling from an actual instrument. 

- Reduced spatial sampling introduces uncertainty into 
seasonal-regional mean AOT assessments. This was 
found for the along-track and across-track MODIS, 
sampling, and again as well for the MERRAero sam- 
pling. Seasonal peaks in the A AOT spatial sampling ar- 
tifact were as large as 0.1 in high AOT regions, repre- 
senting uncertainties as large as 60 % of the full-swath 
AOT values in some places. This uncertainty was some- 
what reduced, but not eliminated, when narrow-swath 
sampling was considered (i.e., when along-track curtain 
sampling was discarded). 

- Reduced spatial sampling introduces uncertainty into 
seasonal-regional mean aerosol radiative forcing as- 
sessments. This analysis was performed in context of 
the MERRAero results, and we found that spatial sam- 
pling - convolving variability in aerosol loading, sur- 
face conditions, and cloud cover - introduced uncer- 
tainty of more than 1 W m“^ and as much as 2-3 W m“^ 
in most regions considered. 

- Spatial sampling affects the derived magnitude and as- 
signment of statistical significance in aerosol trends. 
Along-track curtain sampling results in reduced trend 
magnitude and essentially eliminates statistical confi- 
dence in the derived, decadal-scale trends when the data 
are aggregated at high spatial resolutions. Trend mag- 
nitudes and statistical significance were more similar to 
full-swath values for the narrow-swath sampling. 

- Aerosol trends and statistical significance were found 
to be similar across sampling strategies when the trends 
were composed from coarsely gridded aggregates of the 
sub-sampled MODIS AOT data, suggesting that sin- 
gle pixel width sampling may be sufficient to detect 
and attribute trends at spatial scales on the order of 
1000 km. However, to determine radiative forcing, the 
surface albedo must also be taken into consideration, 
and on these large spatial scales, surface properties, par- 
ticularly over land where AOD tends to be highest, vary 
over much shorter distances. 

In summary, looking at both MODIS data and results of a 
global aerosol reanalysis, we conclude that spatial sampling 
matters. We further conclude that our analysis puts a lower 
bound on the uncertainty in the seasonal-regional AOT intro- 
duced by spatial sampling. Consideration of additional sam- 
plings (e.g., adding more along-track curtains) could act only 
in the direction of increasing the magnitude of our computed 
AAOT sampling artifact. There is also subsequent uncer- 
tainty attached to estimates of DARF due to spatial sampling, 
dependent not only on the variability in the aerosol load but 
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also on surface reflectivity and clouds, and we were able to 
place a lower bound on this through our aerosol model anal- 
ysis. 

The Supplement related to this article is available online 
at doi:10.5194/amt-7-2313-2014-supplement. 
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